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Abstract: We analyze the Snell envelope with path dependent multiplicative 
optimality criteria. Especially for this case, we propose a variation of the Snell 
envelope backward recursion which allows to extend some classical approxima- 
tion schemes to the multiplicatively path dependent case. In this framework, 
we propose an importance sampling particle approximation scheme based on a 
specific change of measure, designed to concentrate the computational effort in 
regions pointed out by the criteria. This new algorithm is theoritically studied. 
We provide non asymptotic convergence estimates and prove that the resulting 
estimator is high biased. 
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Enveloppe de Snell avec des criteres d'optimalite 
multiplicativement dependants de chemin 

Resume : Nous analysons l'enveloppe de Snell avec des criteres d'optimalite 
multiplicativement dependants de chemin. Surtout pour ce cas, nous pro- 
posons une variation de backward recurrence de l'enveloppe de Snell qui permet 
d'etendre certains schemas d'approximation classique a ce cas special. Dans ce 
cadre, nous proposons un schema d'approximation particule d'echantillonnage 
importance base sur un changement de mesure specifique, destine a concentrer 
l'effort de calcul dans les regions soulignees par les criteres. Ce nouvel algorithme 
est theoriquement etudie. Nous fournissons des estimations non convergence as- 
ymptotique et de prouver que l'estimateur resultant est surestime. 

Mots-cles : enveloppe de Snell, option americain, modele particule, evene- 
ments rares 
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1 Introduction 

The Snell envelope is related to the calculation of the optimal stopping time 
of a random process based on a given optimality criteria. In this paper, we 
are interested in some complicated optimality criteria, especially the multiplica- 
tively path dependent case. In other words, given a random process (^ffc)o<fc<n 
and some gain functions (fk)o<k<n and (G k )o< k < n , we want to maximize the 
expected gain E(/ T (X T ) ni=o G k {X k )) by choosing r on a set of random stop- 
ping times T . For example, in finance, the multiplicative optimality criteria 
{Gk)o<k<n El can be interpreted as a discount factor related to a stochastic 
interest rate (taking then an exponential form), or as an obstacle for exotic op- 
tions such as barriers in knock out options (taking then the form of indicator 
functions). 

In the discrete time setting, these problems associated with Snell envelope 
are defined in terms of a given Markov process (X k ) k >o taking values in some 
sequence of measurable state spaces (E n ,£ k ) k >o adapted to the natural filtration 
J- = (J-fe)fe>o- We let 770 = Law(Xo) be the initial distribution on Eg, and we 
denote by M k {xk-i, dxk) the elementary Markov transition of the chain from 
Ek-i into E k . For a given time horizon n and any fc € {0, ...,n}, we let 
Tit be the set of all stopping times r taking values in {k, . . . , n}. For a given 
sequence of non negative measurable functions f k on Ek, we define a target 
process Z k = f k (X k ). Then (U k ) <k<n the Snell envelope of process (Z k ) Q < k < n 
is defined by a recursive formula: 

U k = Z k VE(U k+1 \F k ) 

with terminal condition U„ = Z n . The main property of the Snell envelope 
defined as above is 

U k = sup E(Z T \F k ) = ¥.{Z T .\F k ) with r k = min{fc < j < n : Uj = ZA e T k 

TGTfc k 

Then the computation of the Snell envelope (U k )o< k < n amounts to solve the 
following backward functional equation^] 

u k = fk V M k+ i(uk+i) (1.1) 

for any < k < n with the terminal condition u„ = ,f n . 

But at this level of generality, we can hardly have a closed solution of the 
function u k . In this context, lots of numerical approximation schemes have 
been proposed. Most of them amounts to replace in recursion (jl.ip the pair of 
functions and Markov transitions (/&, M k )o< k < n by some approximation model 
(fk, M k )o< k <n on some possibly reduced measurable subsets E k C E k . In pa- 
per [5], the authors provided a general robustness lemma to estimate the error 
related to the resulting approximation u k of the Snell eveloppe u k , for several 
types of approximation models (/&, Mh)o<k<n- 



In present paper, if not specified, when one talks about the potential rare event Gj. or 
optimality criteria G^, it means G k (Xk) 

2 Consult the last paragraph of this section for a statement of the notation used in this 
article. 
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Lemma 1.1 For any < k < n, on the state space E k , we have that 

n n—1 

|«fc-u*| <^M k4 \f t - f t \ + ^M k:l \{Ai l+1 - M l+1 )u l+1 \ • 

l=k l=k 

This lemma provides a natural way to compare and combine different approxi- 
mation models. In the present paper, this Lemma will be applied in the specific 
framework of a multiplicative optimality criteria. 

Let us come back now to the multiplicatively path dependent case that we 
mentioned in the beginning of the article. Instead of ¥,(f T (X T )) we want to 
maximize E(/ T (X T ) Hp=o G P (X P )) on ^ ne stopping times set T, In this situa- 
tion, a natural way is to consider the path (Xq . . . X k )o< k <n as a new Markov 
chain (X k )o<k<n on path spaces and associate with transitions given for any 
Xk-i = (x ,...,x k -i) 6 {E x- ■ -xE k -i) andxfc = (x' ,...,x' k ) G (E x---xE k ) 
by the following formula 

M k (xk-i,dx'k) = <y x »-i(dXfc-i) M kWk-i,dx' k ) . 

Then, let us denote by Uk(xo ■ ■ ■ x k ) , the Snell envelope defined with a path 
version of recursion (|l.f [) : 

fe-l 

u k (x ,...,x k ) = [fk{x k ) Y[ G p( x p)] V M k+1 (u k+1 )(x , ■ ■ ■ ,x k ) , (1.2) 

p=0 

for < k < n with terminal value u„(xo, . . . ,x n ) = f n (x n ) Tlp^o G p (x p ). At 
this stage, two difficulties may arise. First, the above recursion (|1.2p seem to 
require the approximation of a k + 1 dimensional function at each time step 
from k = n — 1 up to k = 0. Second, when the optimality criteria G p is 
localized in a specific region of E p , for each p, then the product Hp=o G p (x p ) 
can be interpreted as a rare event. Hence, at first glance, the computation of 
Snell envelopes in the multiplicatively path dependent case seem to combine 
two additional numerical difficulties w.r.t. to the standard case, related to the 
computation of conditional expectations in a both high dimensional and rare 
event situation. 

These issues are considered in Section 2, of the present paper. The dimen- 
sionality problem is easily bypassed by considering an intermediate standard 
Snell envelope (v k )o< k < n , without path dependent criteria, which is directly 
related to the multiplicatively path dependent Snell envelope, by the relation 
Ufc(xo, x\...x k ) = rip=o G p (x p )v k (x k ), for all < k < n. Hence, computing the 
original Snell envelope u k can be done by using one of the many approximation 
schemes developped for the standard (non path dependent) case. Then, to deal 
with the rare event problem, we propose a change of measure which allows to 
concentrate the computational effort in the regions of interest w.r.t. the criteria 

n — 1 • 

In Section 3, we propose a Monte Carlo algorithm to compute the multiplica- 
tively path dependent Snell envelope, on the base of this intermediate standard 
Snell envelope under a new equivalent measure defined in the previous section. 
This new approximation scheme is based on the stochastic mesh method intro- 
duced by M. Broadie and P. Glasserman in their seminal paper [3] (see also [7], 
for some recent refinements). The principal idea of original Broadie-Glasserman 
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model is to make a change of probability, under the assumption that the Markov 
transitions Mk(x, ■) are absolutely continuous w.r.t. some other measure rju on 
Eh, with positive Radon Nikodym derivatives Rk(x,y) = (%/)■ But in 

most cases, we do not know the density function of some good choice of r]k- 
So in UJ], the authors provide a variation of Broadie-Glasserman model that 
replaces not only r]k but also the Radon-Nikodym derivatives Rk with the ap- 
proximation model (rff.,Rf-). The model introduced in the present article is an 
extension to multiplicatively path dependent functions. 

In Section 4, the proposed Monte carlo algorithm is theorically analysed us- 
ing an interacting particle system interpretation. We provide non asymptotic 
convergence estimates and prove that the resulting estimator is high biased. 

For the convenience of the reader, we end this introduction with some no- 
tation used in the present article. We denote respectively by V(E), and B(E), 
the set of all probability measures on some measurable space (E,£), and the 
Banach space of all bounded and measurable functions / equipped with the 
uniform norm ||/||. We let fi(f) = J n(dx) f(x), be the Lebesgue integral of a 
function / £ B(E), w.r.t. a measure fj, g ViE). 

We recall that a bounded integral kernel M(x, dy) from a measurable space 
(E,£) into an auxiliary measurable space (E',£') is an operator / H> M(f) 
from B(E') into B(E) such that the functions 

x^M(f)(x) := f M(x,dy)f(y) 

JE' 

are ^-measurable and bounded, for any / S B(E'). In the above displayed 
formulae, dy stands for an infinitesimal neighborhood of a point y in E'. Some- 
times, for indicator functions f = 1a, with A € £, we also use the notation 
M(x, A) := M(1a)(x). The kernel M also generates a dual operator p, n- fj,M 
from M{E) into M(E') defined by ((iM)(f) := n(M{f)). A Markov kernel 
is a positive and bounded integral operator M with M(l) = 1. Given a pair 
of bounded integral operators (Mi,M-2), we let (M1M2) be the composition 
operator defined by (M\M2)(f) = Mi(M2(/)). Given a sequence of bounded 
integral operators M n from some state space E n _i into another E n , we set 
Mk,i '■= Mk+iMk+2 ■ ■ ■ Mi, for any k < I, with the convention Mk,k = Id, the 
identity operator. In the context of finite state spaces, these integral opera- 
tions coincide with the traditional matrix operations on multidimensional state 
spaces. 

We also assume that the reference Markov chain X n with initial distribution 
rjo G V(Eq), and elementary transitions M n (x n -\, dx n ) from E n -\ into E n is 
defined on some filtered probability space (f2, T , P^), and we use the notation 
Ep ?Jo to denote the expectations w.r.t. P, ?0 . In this notation, for all n > 1 and 
for any /„ G B(E n ), we have that 

E P TO {fni.Xn^Fn-i} = M n / n (X„_l) := / M n (X„_l, dx n ) fn(x n ) 

with the a- field T n = (j(Xq, . . . , X„) generated by the sequence of random 
variables X p , from the origin p = up to the time p — n. We also use the 
conventions JTg = 1, and = 0. 
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2 Snell envelope with multiplicatively path de- 
pendent functions and change of measure 

Suppose (Xk)o<k<n is a Markov chain on continuous state spaces (E k , £k)o<k<n 
with an initial distribution r/o on Eq, a collection of Markov transitions Mk(xk-i, 
dxk) from E k _\ to Ek and a given final time horizon n. We also assume that the 
chain X k is defined on a filtered probability space (ft, J 7 , P). In this situation, 
the historical process X k := (Xo, . . . ,X k ) can be seen as a Markov chain with 
transitions given for any Xk-i = {xq, . . . , x k -i) G Eq x • • • x £7fc_i and xj,, = 
(x' Q , . . . , x' k ) £ Eq x • • • x Ek by the following formula 

M k (Xk-i,dx'k) = ^xk-iWk-i) M k (x' k _ x ,dx' k ) . 

We denote by (Pfc)o<fc<n a sequence of probabilities of path (X k )o< k < n . For a 
given collection of real valued functions (f k )o<k<n and (Gk)o<k<n, defined on 
(Ek)o<k<m we define a class of real valued functions (F k )o< k <n defined on the 
product spaces (E x ■ ■ ■ x E k ) < k < n by 

F k (x ,- ■■ ,Xk) ■■= fk(xk) G p (x p ) , for all < k < n . 

0<p<fc-l 

To maximize the expected gain E(F T (X T )) w.r.t. r in a set of random stopping 
times T, one is interested in computing the Snell envelope (ufc)o<fc<n associated 
to the gain functions (F k )o <k < n and solution to the following recursion 

{^n (xq , , X n ^ — F n (xQ, , X n *) 
u k (x , ■ ■ ■ ,x k ) = F k (x , ■■■ ,x k )V M k+ i(xi k+ i)(x , . ..,x fe ),V < k < n - 1 

(2.1) 

Now, let us consider the standard (non path dependent) Snell envelope (v k )o< k < n 
associated to the gain functions ( f k )o< k < n and satisfying the following recursion 

{^n{Xn) — fni^Xji) 
v k (x k ) = f k (x k ) V [Gk(xk)M k+1 (vk + x)(x k )] , for all < k < n - 1 . 

(2.2) 

For all < k < n, let us denote by the real valued functions defined on 
E x---xE k , such that v k (x , ■ ■ ■ ,x k ) := v k (x k ) Ilp=d G p ( x p ) ■ B ^ construction, 
one can easily check that for all < k < n, u k = v k and in particular Uo(a;o) = 
vq(xq). Indeed, one can verify that (~v k )o<k<n follows the same recursion (|2.1[) 
as (ufe)o<fe<n- First, we note that they share the same terminal condition, 

n — 1 n— 1 

v„(x , ■ ■ ■ ,x n ) = v n (x n ) Y[ Gp(x p ) = f n (x k ) G p (x p ) 

Fn (xq , • • • , X n J — \l n (xq , • • • , X n ) . 
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Then at time step k, we observe that they follow the same recursion 



v fc (x , • • • ,x k ) 
fe-i 

Vk{xk) Y\_ Gp{x p ) 
fc-i 

fk{xk) Y[ g p( x p) 

p=0 



fc-1 



M k+1 (xk,dxk+i)vk+i{x k+ i)Gk(x k ) Y\_ Gp(x p ) 

p=0 



= F k (x , ■ ■ ■ ,x k ) V A^fc + i(v fc+ i)(a;o, . . . ,x k ) ■ 

Now that we have underlined the link between M k and v k , our aim is then 
to compute the latter. The recursion (|2.2p implies that it is not relevant to 
compute precisely the conditional expectation Mk+\(v k +i){x k ) when the value 
of the criteria Gk{xk) is zero or very small. Similarly, notice that Vk+i is likely 
to reach high values when Gk+i does, hence from a variance reduction point 
of view, when approximating the conditional expectation Mk+\{vk+i){xk) by a 
Monte Carlo method, it seems relevant to concentrate the simulations in the 
regions of E k +\ where Gk+i reaches high values. Hence, to avoid the potential 
rare events G, we propose to consider the following change of measure one the 
measurable product space (Eq x • • • x E n , £q x • • • x £„), 



1 



T[G k 

k=0 



with Z n =E[Y[G k (X k ) )=]Jvk(G k ), 



\k=0 



fc=0 



(2.3) 

where r\k is the probability measure defined on Ek such that, for any measurable 
function / on E k 

^(f(Xk)Up^ G p (X p )) 



Vk{f) ■-- 



The measures (?7fc)o<fc<n defined above can be seen as the laws of (^fc)o<fc<n 
under probability Q. Loosely speaking, the process (Xfc)o<fc< n with distribution 
(i]k)o<k<n is designed under the constrain {Yl p= o Gp)o<k<n- An intuitive inter- 
pretation comes by setting Gk{xk) = lA k {xk) with A k C Ek, then the process 
with distribution rjk is just the ones surviving in the subsets Ak- It follows that 
the measures rjk seem to be a relevant choice for the change of probability in 
our path dependent situation. 

Furthermore, it is also important to observe that, for any measurable function 
/ on E k 

Vk(f) = tt; ^ ■ (2.4) 

T)k-l(Gk-l) 

We denote the recursive relation between rjk and 7]k-i by introducing the oper- 
ators $fe such that, for all 1 < k < n 

Vk = $k(%-i) • (2.5) 
Let us now introduce the integral operator Qk such that, for all 1 < k < n 



h(f)(x k -i) 



Gk-i(xk-i)M k {x k -i,dxk)f{xk) 



(2.6) 
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In further developments of this article, we suppose that M k (x k —i, ■) are equiva- 
lent to some measures X k , for any < k < n and x k -i G i.e. there exists 
a collection of positive functions H k and measures such that: 

M k (x k -i,dxk) = H k (x k -i,Xk)\k(dx k ) . (2.7) 

Now, we are in a position to state the following Lemma. 

Lemma 2.1 For any measure r\ on E k , recursion h2.2\) defining v k can be 
rewritten: 

Vk(x k ) = fk{xk) V Qk+i{vk+i){xk) = fk(xk) V $fe+i(7y) r^ +l( ?; ^ fc+i 



V d<Z> k +i(v) 



for any Xk £ Ek, where 



dQ k +i(x k , ■) . _ G k {x k )H k+ i (x k , x k +i )y(Gk) 
d$ k+1 ( V ) [Xk+l) ~ r)(G k H k+1 (.,x k+1 )) 

for any (x k ,x k +i) € E k x E k+X . 
Proof: 

Under Assumption (|2.7[) . we have immediately the following formula 

M k+1 (x k ,dx k+ i) = H k+1 (x k ,x k+ i) — — ^ k ( G k) —n k+ i{dx k+ i) . (2.8) 

r] k (G k H k+ i(-,Xk+i)) 

Now, note that the above equation is still valid for any measure rj, 

f](G ) 

M k+1 (x k ,dx k +i) = Hk+i(x k ,Xk+i) 7 —$k+i(r])(dxk+i) ■ 

r]{G k H k+1 (-,Xk+i)) 

(2.9) 

Hence, the Radon Nikodym derivative of M k +i(xk, dxk+i) w.r.t. <& k+ i{j]) is 
such that 

dM k+1 (x k ,-) rj{G k ) 

— — — — (Xfc+l) = Hk+l(x k ,Xk + l) rr . 2.10 

d^k+i{V) n(G k H k+1 (-,x k+1 )) 
We end the proof by applying above arguments to recursion (|2.2D . 



3 A particle approximation scheme 

From the above discussion, we conclude that the distributions (r) k )o< k < n are a 
very good choice for the change of probability for the stochastic mesh model. In 
this section, we first propose a particle model to sample the random variables 
according to these distributions, then we describe the resulting particle scheme 
proposed to approximate the Snell envelope (v k )o<k<n- 
By definition (|2.5[) of we have the following formula 

®k(r)k-i) = i]k-iK k ^ lk _ 1 = <q k -xS k -x !Vh _ l M k = y Gk _ 1 (rik-i)M k . (3.1) 
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Where K k>Vk _ 1 , S k -i^ kl and ^c k _i are defined as follows: 
K k , Vk _ 1 (x k -i,dx k ) = (S k - 1>Vk _ 1 M k )(x k ^ 1 ,dx k ) 

= J S k -i, Vk _ l {xk-i,dx' k _ l )M k (x k _ l ,dx k ) , 
(x, dx') = eG k - 1 (x)S x (dx') + (1 - eG k - 1 (x))'f Gk _ 1 (n^idx 1 ) 

VG^AVk-iXdx) = ^-^ m-iidx) , 

where the real e is such that eG takes its values [0, 1]. 

More generally, the operations \1/ and S can be expressed as ( 7 ?)(/) = ifoT) ~ 
v S v (f) with £„(/) = eGf + (1 - cG)*g(»7)(/). 

The particle approximation provided in the present paper is defined in terms of 
a Markov chain £^ = (£ k )i<i<N on the product state spaces , where the 
given integer N is the number of particles sampled in every instant. The initial 
particle system, (jf^ = ^£o*'^) , is a collection of N i.i.d. random copies 

of Xq. We let J 7 ^ be the sigma- field generated by the particle approximation 
model from the origin, up to time k. To simplify the presentation, when there 
is no confusion we suppress the population size parameter N, and we write £ k 
and £|. instead of £ k N ^ and £ k ' N ^ . By construction, £ k is a particle model with 
a selection transition and a mutation type exploration i.e. the evolution from 
£,k to ^fe+i is composed by two steps: 

, r Selection ^ \ » r Mutation » r 

& G < ► & := (a) „ _ e £f ► 6 + i G Sf+i • (3.2) 



k.rj 



Then we define r/^ and 77^ as the occupation measures after the mutation and 
the selection steps. More precisely, 

E % «-* €--=± E 

l<j<AT l<i<jV 

During the selection transition S k r/ N, for < i < N with a probability eG k (£ k ) 

we decide to skip the selection step i.e. we leave stay on particle and with 
probability 1 — eG k (£ k ) we decide to do the following selection: Q. randomly 

takes the value in fi for < j < N with distribution ^w^^tti ; ■ Note that 

when eGfc = 1, the selection is skipped ( i.e. = £ k ) so that the model 
corresponds exactly to the Broadie-Glasserman type model analysed by P. Del 
Moral and P. Hu et al. [5]. Hence, the factor e can be interpreted as a level of 
selection against the rare events. 

During the mutation transition ~-+ £ k +i, every selected individual evolves 
randomly to a new individual C k +i = £ randomly chosen with the distribution 
M k+1 (C k ,dx), for 1 < i < N. 

It is important to observe that by construction, rj k+1 is the empirical measure 
associated with N conditionally independent and identically distributed random 
individual C k +i with common distribution 1(77^). 
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Now, we are in a position to describe precisely the new approximation scheme 
proposed to estimate the Snell envelope (vk)o<k<n- The main idea consists in 
taking rj = 77^, in Lemma 12. 11 then observing that Snell envelope (vk)o<k<n is 
solution of the following recursion, for all < k < n, 

Now, if $>k+i (Vk) i s wen estimated by it is relevant to approximate Vk by 

Vk defined by the following backward recursion 

Vn fn 



v k (x k ) = fk{xk) V nj^ +1 ( ^Tw gHi) for all < fc < n , 

(3.3) 

Note that in the above fomula (|3.3[) . the function Vfc is defined not only on 
but on the whole state space E k . 
To simplify notations, we set 

d<S>k+i{Vk) 

Finally, with this notation, the real Snell envelope (vk)o<k<n and the approxi- 
mation (ik)o<k<n are such that, for all < k < n, 

Vk = fk V Q k+1 (v k +i) 
Vk = fk V Qk+i(vk+i) ■ 



4 Convergence and bias analysis 

By the previous construction, we can approximate <f>k+i(Vk) by Vk+i- ^ n 
this section, we will first analyze the error associated with that approximation 
and then derive an error bound for the resulting Snell envelope approximation 
scheme. To simplify notations, in further development, we consider the random 
fields V k N defined as 



The following lemma shows the conditional unbiasedness property and mean 
error estimates for the approximation n k+1 of 

Lemma 4.1 For any integer p > 1, we denote by p' the smallest even integer 
greater than p. In this notation, for any < k < n and any integrable function 
f on Ek+i, we have 

E(% W +1 (/)|-Ff)=<i> fe+1 ( % Ar )(/) 

and 



E 



(\V k N {f)\ P \^)' < 2 a(p) [* k+l (v?W\ p ') 



with the collection of constants 

a(2p) 2 P = (2p)_ 2-p and a(2p + l) 2 ^ 1 = (2 p+1)p+ 1 2-^+ 1 ^ 

y/p + 1/2 
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Proof : The conditional unbiasedness property is easily proved as follows 

1 N 



N 

i=l 



1 N 

= (r ? fif fc+1 , < )(/) = ^ +1 (^)(/). 
Then the above equality implies 

E (I - W^)] (/)f 1^) * < E (I fo&i - (/)| P l<) * . 

where /xj^ := -k J^ i=1 <5y^ stands for an independent copy of given 77^. 
Using Khintchine's type inequalities yields that 

jN*(\[r&i-l&M'\r?)* < 2a( P )E(\f(e k+1 )f\^f 

= 2a(p) [$ fc+1 (r ? f)(|/r')]^ . 
We end the proof by combining the above two inequalities. 

■ 

A consequence of the unbiasedness property proved in Lemma 14.11 is that 

E(Q k+1 (f)(x k M) = Q k+1 (f)(x k ) . 

To estimate the error between v k and the approximation v k , it is usefull to intro- 
duce the following random integral operator R k such that for any measurable 
function on E k+ i, 



«2ui(/)(a*) - Vn (Q k+ x(f)(x k ) - Qjt+i(/)(x fc )) 



Note that 



then, applying again Lemma 14 . 1 1 implies the following Khintchine's type inequal- 



ity 

n\R k+ i(v 



< 2 a(p) 



®k+i(Vk)( dx k+i) ( ^t+i^O ^ Xk+l ^ fc+1 ^ Xk+1 ^ 



Let Q k j = Qfc+iQfc+2 • ■ • Qi for any < k < I < n, then it follows easily, by 
recursion, that 

^(Qk,i(f)(xk)\Vk) = Qk,i(f)(xk) ■ 
Now, by Lemma 1 1.11 we conclude 

VN\(v k -d k )\< QmI«i)(^+i)I • (4.1) 

fc<Z<ri 

We are now in position to state the main result of this paper. 
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Theorem 4.2 For any < k < n and any integer p > 1, we /iai>e 



sup ||(55]fe — Ufc)(aj)||i < 



2 o(p) 



k<l<n 



■Qk,l 



}k,i+i(hf +1 1 vf +1 )(x) 



with a collection of constants q^ i and functions h k defined as 

p'-i 



qk,i 



i-i 

\Gi\\ \\h k+ i\\ J] II G * 

m—k 



and h k (x k ) 



sup 



Hk(x,Xk) 



Proof : First, decomposition (|4.ip yields 

ViV||(^- Wfc )(x)|j ip < ^ |Q fe ,;|« 1 )(^+i)|(^ 



k<l<n 



x,ye-E fc _i H k (y,x k ) 

(4.2) 



for all ie G -Ei- 



Note that 



i-i 



IIQfc,z(l)ll < & fc.f , where 6 fe . ( := ||/ifc+i|| ||G m || . 

m—k 

Then it follows easily that for any integrable function / on Ei 

(Qk.i(f)) p < {h,i) p ~ 1 Qk,i(f p ) ■ 

This yields that 

QkJ |(-Rt+i))(«I+i)| (x) <(6 M )^ i E(QM(|«i))(t' i+ i)|) P (^) 

Lin \ 



Applying Lemma 14.11 to the right-hand side of the above inequality, we obtain 
for any xi £ E; 



< 2 o(p) 



Qi+iinF ) (d^+i ) ( ^v? ( a; '+i (^+0 



from which we find that 



(0(«w)(«i)rw 



< 2 o(p) 



Qi + i(a;;,da;/ + i) 



/ dQi +1 (xi , •) 



Am 



By definition (|4.2[) of functions hi+i and in developing the Radon Nikodym 
derivative, we obtain 



dQi+i{xi,-) 



T)( f (Gi)G l {x l )H l+1 (x l ,xi 
r,f(G l H l+1 )(-,x l+1 ) 



-fl; 



< ||G/||/iz + i(xi + i) , 
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which implies 

E(i« 1 ))(« l+ i)(x J )ri^)' 

<2a(p)\\G l \\^ , 



x{xudxi +1 ) (hi +1 (x l+1 )) v 1 vi +1 {xi + i) p 



Gathering the above arguments, we conclude that 

\\(v k -V k )(x)\\ Lp < ^ 2 %r 1k,l (Qk.l+lihf+^V^^y 



k<l<n 



Remarks : The constants qk i could be largely reduced. In fact, qk i comes from 
bounding || n m rtt(G m )\\ Lp . In g], the authors proved || ]J m G m \\ L ' 2 + ^f^ as 
a non asymptotic boundary for |j FT T]^ t (G m )\\L 2 - In most cases, the functions 
G take their values in [0,1], then the majoration linm^mll — 1 holds, but 
II Ilm G m \\ L2 is very small. 

When the function G vanishes in some regions of the state space, we also 
mention that the particle model is only defined up to the first time t n = k such 
that rj^ (Gfe) = 0. We can prove that the event {t n < n} has an exponentially 
small probability to occur, with the number of particles N. In fact, the estimates 
presented in the above theorems can be extended to this singular situation by 
replacing Vk by the particle estimates Vkl T N >n . The stochastic analysis of these 
singular models are quite technical, for further details we refer the reader to 
section 7.2.2 and section 7.4 in the book [S]. 

To understand better the L p -mean error bounds in the above theorem, we 
deduce the following exponential concentration inequality 

Proposition 4.3 For any < k < n any and any e > 0, we have 



sup P 



(\v k (x) - v k {x)\ >^=+e^j< exp (-Ne 2 /c 2 )) , (4.3) 



with constant c = J2k<i<n 2 1k,i (Q/u+itMVi^'f+iH 2 '')) " ■ 

Proof : This result is a direct consequence from the fact that for any non 
negative random variable U such that 

3b < oo s.t. Vr > 1 E {U r Y < a(r) b => P (U > b + e) < exp (-e 2 /(2b 2 )) . 
To check this claim, we develop the exponential and verify that 

Vt > E (e fU ) < exp ( + bt] ^P(U>b + e)< exp (- sup (et - 



Simarly to the orginal Broadie-Glasserman model, the following proposition 
shows that in this model we also over-estimate the Snell envelope. 
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Proposition 4.4 For any < k < n and any x k G E k 

E(v k {x k )) > v k (x k ) . (4.4) 

Proof: 

We can easily prove this inequality with a simple backward induction. The ter- 
minal condition v n = v n implies directly the inequality at instant n. Assuming 
the inequality at time k + 1, then the Jensen's inequality implies 



E(v k {xk)) > fk(xk)VE[Q k+1 Vk +1 (x k ) 

= / fc (a:fc)VEj / Q k+1 (x k ,dx k+1 )E(G k+1 (x k+1 )\F? +1 ) 
By the induction assumption at time k + 1, we have 

E / Q k+1 {x k ,dx k+ i)E(v k+ i{x k+ i)\Fk +1 )\ > E (Q k+1 v k+1 (x k ) 

= Qk+iv k+ i(x k ) . 
Then the inequatily still holds at time k, which completes the proof. 
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